Field induced ordering in highly frustrated antiferromagnets 
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We predict that an external field can induce a spin order in highly frustrated classical Heisenberg 
magnets. We find analytically stabilization of collinear states by thermal fluctuations at a one- 
third of the saturation field for kagome and garnet lattices and at a half of the saturation field for 
pyrochlore and frustrated square lattices. This effect is studied numerically for the frustrated square- 
lattice antiferromagnet by Monte Carlo simulations for classical spins and by exact diagonalization 
for S — 1/2. The field induced collinear states have a spin gap and produce magnetization plateaus. 

PACS numbers: 75.10.Hk, 75.10.Jm, 75.50Ee, 75.40.Mg 



Frustrated magnets, classical and quantum, exhibit 
spectacular and often unexpected behaviors [Q . A few of 
them do not order at any T > 0. Extensive residual en- 
tropy in the classical ground state in such cases prevents 
the realization of the usual order by disorder scenario Q . 
An external magnetic field changes the degeneracy and 
topology of the ground state manifold stabilizing, even- 
tually, the (nondegenerate) saturated state at H > H sat - 
If the order by disorder effect occurs in a finite field and 
suppresses residual entropy, then, such an effect can be 
used for practical applications. During a demagnetiza- 
tion process a spin system has to regain its entropy and, 
therefore, the whole crystal will cool down. Examples of 
geometrically frustrated AFMs on kagome |j| , pyrochlore 
and garnet || lattices include magnetic compounds 
with rather small exchange constants opening the way 
for experimental tests of their finite field behavior. 

In this Letter, we predict that fluctuations stabilize 
collinear spin configurations at rational values of the field 
H/H sat = 1/2 or 1/3. We present an analytical proof of 
the field induced ordering driven by thermal fluctuations 
and suggest a similar role for quantum fluctuations on 
the basis of spin- wave and numerical results. In fact, an 
ordered spin phase in a finite field has already been ob- 
served in gadolinium garnet GdaGasO]^, though dipolar 
anisotropy plays a crucial role in this material |^| . 

To be specific we use as an example the frustrated 
square-lattice antiferromagnet (FSAFM) ||] in an exter- 
nal field: 

n = j£s,-s i + j' ]T s, • s ; // £ s; . (i) 

n.n. n.n. n. i 

First, we present analytical arguments, making them as 
general as possible, in order to include the other highly 
frustrated AFMs. Second, we consider numerical results 
for classical and quantum FSAFMs in a magnetic field, 
which confirm our predictions. 

Let us briefly discuss magnetically ordered phases of 
the FSAFM. For small diagonal exchange J' < 0.5 J, clas- 
sical spins form the Neel state in zero field. At J' > 0.5 J 



the classical ground state consists of two \/2 x \/2 inter- 
penetrating antiferromagnets, which are locked by fluc- 
tuations in a striped AFM state described by a single 
wave- vector (n, 0) or (0, n) An applied magnetic field 
cants the spins and creates an easy plane for the AFM 
sublattices. Thus, both the Neel and striped states have 
a transverse spin order in a finite field. Thermal and 
quantum fluctuations destroy the transverse order in the 
vicinity of the highly frustrated point J' = 0.5J. As 
an illustration we present in Fig. ^ the phase diagram 
for the spin-1/2 model at T = obtained in the linear 
spin-wave theory. This phase diagram agrees with the 
previous zero-field studies B, which suggest a quantum 
disordered ground state for S = \ near J' = 0.5J. Since 
the disordered singlet phase has a spin gap, it becomes 
unstable above a finite critical field. 
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FIG. 1. Zero temperature phase diagram of the spin-1/2 
FSAFM in external field obtained from the linear spin-wave 
theory. The shaded region indicates the stability region of the 
quantum disordered state with a spin gap. 

We now focus on the classical critical point J' — 0.5 J, 
where the Hamiltonian (|l|) can be written up to a con- 
stant term as a sum over edge-sharing plaquettes: 

^ = F^{ J|LQ|2_jff ^}' (2) 
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where L Q = J^iea &i 1S t ne total spin in plaquette a and 
n = 2. The block form of the spin Hamiltonian (||) is 
common to all highly frustrated spin models. For AFMs 
on kagome and garnet lattices the blocks are triangles 
and for a pyrochlore AFM they are tetrahedra with the 
corner-sharing arrangements and n = 1 in all three cases. 

A spin configuration minimizes the energy (|^) at H = 
provided L Q = for each plaquette. This constraint can 
be satisfied for many classical states. We estimate the 
number of continuous degrees of freedom in the ground 
state of the A-site FSAFM as ~ A 1 / 2 . For pyrochlore 
and kagome AFMs the degeneracy is larger and the di- 
mensionality of the ground state manifold scales with A 
P,H0]. Thus, zero-field properties are not universal for 
disordered frustrated AFMs. We now show that a uni- 
versal behavior does appear in a magnetic field. 




FIG. 2. The collinear spin structure for a four-spin block 
at H = 0.5-ff sa t. Rotations by angles 9 and tp parameterize 
the zero-energy and soft modes in this state. 

In a finite field the classical energy (Q) is minimized 
for spin configurations with the plaquette magnetiza- 
tion L z a = H/(2J). There are many degenerate clas- 
sical states which satisfy this constraint below the sat- 
uration field -ff sa t = 8JS (pyrochlore and FSAFM) or 
6JS (kagome and garnet). Generally, all these states are 
noncollinear as, e.g., canted Neel and striped phases for 
FSAFM. Collinear spin arrangements appear only at spe- 
cial rational values of the applied field: at H c — ^H sat for 
4-spin blocks (FSAFM, pyrochlore), the up- up- up-down 
(uuud) structure, see Fig.||[ and at H c = ^H sat for 3-spin 
blocks (kagome, garnet), the up- up-down structure. The 
distribution of down-spins on the lattice in a collinear 
state is not unique: it only obeys the one down-spin per 
plaquette condition. 

Zero-energy modes correspond to continuous distor- 
tions of a given classical ground state, which do not vi- 
olate the magnetization constraint. To construct a zero 
mode for the collinear states described above one has to 
draw an (open) line through the lattice points with adja- 
cent sites on a line occupied by antiparallel spins with no 
plaquettes crossed more than once. A zero mode consists 
of a simultaneous rotation of all spins along the line by 
an angle 9 (Fig. ^). The complete set of zero modes is 
constructed when a set of such lines goes through all up- 
spins. Zero modes connect a uuud state to noncollinear 
ground state configurations. Hence, the collinear states 
are singular points on the ground state manifold. There 
are no other collinear states in the local neighborhood of 



a selected collinear state. 

Moessner and Chalker have recently discussed the 
order by disorder phenomenon in zero field for classi- 
cal frustrated models in terms of a local topology of the 
ground state manifold. They reached definite conclusions 
only for an XY magnet on the pyrochlore lattice and for 
a Hcisenberg kagome lattice AFM. We show that a simi- 
lar analysis predicts stabilization of collinear states near 
special rational values of the applied field. Let x denote 
the coordinates on the ground state manifold and y be 
the transverse directions, which span the rest of the con- 
figuration space. Low-energy excited states are described 
by the quadratic Hamiltonian Ti.2 — X)z £ (( x )?// 2 resulting 
in a probability distribution 



Z(x) - J d yi e-^ cx UlT/efa)] 



1/2 



(3) 



over the ground state manifold at low temperatures. At a 
special point xo some of the stiffnesses e;(xo) may vanish 
making Z(xq) divergent. The corresponding coordinates 
yi describe soft modes. The appearance of the order by 
disorder depends on a number of soft modes, which exist 
for a given classical ground state x . 

A special feature of the collinear states with finite mag- 
netization is that each zero mode generates exactly one 
soft mode. They are constructed in the following way. 
First, all spins along a zero mode line are rotated by 9 
about the axis perpendicular to the field and, second, all 
down spins on the same line are rotated by an angle ip 
about the field direction, see Fig. ||. The total spin of 
each plaquette in a deformed state is L Q w 5(0, 6<p, 2) 
for small 9 and ip and contributes to an energy increase 
~ 9 2 ip 2 . Thus, we identify x\ with 9 and yi with ip. 

To allow the order by disorder selection, not only has 
the probability density Z(xo) to diverge, but the statisti- 
cal weight J Z(x)dx must be concentrated entirely near 
the collinear spin states. We check this by integrating 
Z(x) in the £>-dimensional neighborhood of a collinear 
state parameterized by zero modes. Since each zero mode 
has one soft mode with e s (x) ~ x 2 , the integral 



Z(x)ci> 



(4) 



diverges independently of the actual value of D. There 
is a vanishingly small probability of finding a spin sys- 
tem in one of the noncollinear states surrounding a given 
collinear configuration. Hence, the order by disorder se- 
lection occurs: thermal fluctuations stabilize a discrete 
set of collinear states. This result relies only on the spe- 
cial block structure (||) of the spin Hamiltonian and the 
symmetry properties of the collinear states. Therefore, 
the order by disorder in an external field is a universal ef- 
fect and occurs for all frustrated Heisenberg AFMs: apart 
from the FSAFM the uuud states at H — ^H sat are sta- 
bilized for a pyrochlore AFM and the uud states appear 
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for Heisenberg AFMs on kagome and garnet lattices at 
H = ^H sat - A subsequent selection between collinear 
states with different patterns of down-spins is made to 
maximize D in Eq. (Q) and corresponds in the case of 
FSAFM to the most symmetric q = uuud state. The 
above proof of the field induced ordering is based on the 
special property of the collinear states with finite mag- 
netization. Namely, an equal number of zero and soft 
modes. For a Heisenberg magnet on the pyrochlore lat- 
tice in zero field the same equality holds only approxi- 
mately: in the leading order in the number of spins N. 
Therefore, no analytical conclusion has been reached in 
a zero field case, while numerical simulations indicated 
no ordering {jj. 

For quantum models in the limit S > 1 we check the 
relative stability of different states by comparing their 
zero-point oscillation energies: \ u>k. The magnon 
spectrum of the q = uuud state has four modes: 



wi,2 = JSy/rj^ , w 3j4 = 2JS± JSy/A- 



(5) 



with ?7k = (1 — cosfc x )(l — cosfcj,). Zero-point contribu- 
tions to the energies of the q = uuud state, the Neel 
state, and the striped AFM at J' = 0.5 J and H = ^H sat 
are 0.57JS, 0.69JS", and 0.81J5, respectively. (They 
have the same classical energy.) Thus, this compari- 
son suggests that quantum fluctuations also select the 
collinear states in a magnetic field due to their large num- 
ber of soft modes. Collinear states preserve the 0(2)- 
rotational symmetry. Hence, their renormalized magnon 
spectrum becomes gapped and plateaus arise on the mag- 
netization curve. The spin-density wave of the q = 
uuud state in FSAFM is a superposition of spin harmon- 
ics with wave-vectors (0,0), (7r,7r), (tt, 0) and (0, 7r). 

To test these predictions we have performed Monte 
Carlo (MC) simulations for the model Eq. (]l|) with unit 
vector spins. Fig. [| top panel shows the magnetization 
curve at T = 0.1 J obtained with 2 x 10 5 MC steps per 
spin for each point. The field derivative of the magne- 
tization has a dip around H = 4 J = ^-^sati which in- 
dicates the presence of a new plateau phase. A jump in 
the magnetization and a hysteresis of about AH ~ 0.3J 
suggests a first-order transition to a high-field phase. To 
determine the nature of the spin state at the plateau we 
calculated different components of the static structure 
factor S a ^(q) = pE r!i e iq ' x (S r a Sf +I ). The Neel and 
the collinear states have nonzero transverse components 
gxx _ gyy at q_ ^ ^ anc j q _ ^ 7F) o) or (0, 7r), respec- 
tively. The uuud state has nonzero longitudinal com- 
ponents S zz (q) at all the above vectors simultaneously. 
The static structure factor presented in Fig. ||, bottom 
panel was obtained by averaging 50 'instant shots' sep- 
arated by 10 3 MC steps. In the region of weak (strong) 
diagonal exchange J' the data clearly support the Neel 
(striped) type of spin correlations. Nonzero harmonics 
in the longitudinal structure factor both at q = (n, 0) 



and q = (tt, tt) exist only for 0.494 < J'/ J < 0.508. We 
have also checked that neither a lower field of H = 2 J, 
nor a higher field of H = 6 J induces a nonzero value of 
S zz (q) at these points in the Brillouin zone. Thus, these 
results unambiguously identify the spin configuration on 
the plateau as the q = uuud state. 




FIG. 3. Monte Carlo results for the classical FSAFM: the 
magnetization curve for J'/J = 0.5 (top panel) and the 
static structure factor S a/3 (q) vs. the frustration parameter 
for H = 4J (bottom panel) on a 24x24 lattice at T = 0.1 J. 

The uuud phase breaks the translational symmetry in 
such a way that the wave-vectors belonging to different 
irreducible representations of the space group are mixed: 
(tt,tt) vs. (tt, 0) and (0,7r). Hence, the spin structure of 
the uuud phase is a mixture of different order parameters. 
In such a case the phase transition to a disordered spin 
liquid state at the high-field end of the plateau, Fig. ||, 
must go either via a first order transition or in several 
steps with an intermediate supersolid state. 

We have also studied the quantum spin-1/2 model (|l|) 
at T = by Lanczos diagonalizations of finite clusters. 
Fig. [| presents the magnetization m = M/S vs. field 
at J'/J = 0.6. There are no magnetization plateaus in 
the thermodynamic limit for m > 1/2. To determine 
whether the plateau at m — h remains after finite size 
scaling we show its width in the inset for three cluster 
sizes as a function of J'/J. If a plateau disappears in 
the thermodynamic limit, its width should decrease as 
Am = l/N. Therefore, the ratio of the plateau- widths 
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of the 4x4 and 6x6 clusters should be 16/36 or less if 
the magnetization curve has a nonzero slope for N — > oo. 
The ratio 16/36 is exceeded for m — 1/2 only in the 
interval 0.49 < J'/ J < 0.66, where we would expect to 
have a plateau on an infinite lattice. 




FIG. 4. Exact diagonalization results for the spin-1/2 
FSAFM. Top panel: the magnetization curve for J'/J = 0.6. 
Inset: width of the finite-size plateau with m = | as a func- 
tion of J'/J. The lines are for lattices of size 4x4 (dashed), 
6x4 (dotted), and 6x6 (full). The bold line is an extrapolated 
curve. Bottom panel: static structure factor as a function of 
J'/J atm=iona6x6 cluster. S zz (0, it) is shown by '+', 
S zz (n,n) by 'x', S xx (0, n) by and S x * (tt, tt) by squares. 

Quantum fluctuations have a strong effect on the sta- 
bility of the m = | plateau. Both, the width of the 
plateau for the spin-1/2 system and the range where it 
appears are larger than for the classical model. In addi- 
tion, the parameter range for the plateau is shifted asym- 
metrically around the classical critical point J' = 0.5J. 
The plateau is most pronounced for J' 0.6 J. This is 
the value we used for the presented magnetization curve. 
The bold line in Fig. g was obtained by connecting the 
mid-points of the steps for the largest available systems 
size, except for the plateau at ra = |, for which the cor- 
ners of the 6x6 cluster data were used. The value of the 
spin gap (i.e. the boundary of the m = plateau) was 
taken from Kotov et al. ||. Fig. [|also shows the peaks of 
the static structure factor vs. J'/ J for the m — \ plateau. 
The peaks in S zz (q) , which indicate the presence of the 
uuud state, exist for 0.51 < J'/J < 0.67. 

The discussed mechanism for the magnetization 
plateaus in highly frustrated AFMs has to be contrasted 



with that of weakly coupled ID spin systems (see, e.g., 
S). There, magnetization plateaus correspond to dis- 
ordered states with a gap, that are stable under a weak 
higher-dimensional coupling, while intermediate gapless 
regions are immediately ordered at sufficiently low tem- 
peratures once a higher-dimensional coupling is switched 
on. Here, we found that an external field induces a long- 
range collinear ordering with an excitation gap on the 
plateau, whereas intermediate gapless regions remain dis- 
ordered. Further exploration of such possibilities should 
also be important for the general problem of classifying 
magnetization plateaus in two dimensions p . 
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